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Abstract 

A unified semiclassical time propagator is used to calculate the semiclassical 
time-correlation function in three cartesian dimensions for a particle moving 
in an attractive Coulomb potential. It is demonstrated that under these 
conditions the singularity of the potential does not cause any difficulties and 
the Coulomb interaction can be treated as any other non-singular potential. 
Moreover, by virtue of our three-dimensional calculation, we can explain the 
discrepancies between previous semiclassical and quantum results obtained 
for the one-dimensional radial Coulomb problem. 

3.65.Sq,3.65.G,31.50 



Typeset using REVTgX 



Semiclassical propagation in time has been studied intensively in two dimensions |]-f3- 
There are by far not as many applications to higher dimensional problems, in particular not 
in connection with the singular Coulomb potential. Our motivation for this study is three- 
fold: Firstly to see, if the advanced semiclassical propagation techniques in time, namely 
the Herman-Kluk propagator [§-0], can be implemented for realistic problems of scattering 
theory involving long range forces. Secondly to see, if we can avoid to regularize the Coulomb 
singularity in the classical equations of motion if we work in three (cartesian) dimensions, 
and thirdly, to clarify the reason for the small, but pertinent discrepancies with the quantum 
result in two previous, one-dimensional semiclassical calculations of the hydrogen spectrum 
from the time domain @,^]. As it will turn out, the Coulomb problem with the Hamiltonian 
(we work in atomic units unless stated otherwise) 

#=- + ^ (1) 
2 |r| K J 

can be propagated in time semiclassically without taking any special care of the singularity 
in the potential which poses a lot of difficulties for the one-dimensional radial problem if 
Z < 0, i.e., if the potential is attractive as in the case of hydrogen (Z = —1) which we 
take as an example in the following. The relevant information in the time domain is the 
autocorrelation function 

c(t) = (V>|A#> (2) 

where 

K(r,r',t) = (r\e- im/h \r') (3) 

is the propagator in the coordinate representation. By diagonalizing K in Eq. (Q) one can 
express the autocorrelation function with the time evolution operator U(t), 

C (t) = (mmt)mo)) = (mmt))- w 

This form has the obvious interpretation of correlating the time evolving wavefunction ijj(t) 
at each time with its value at t — 0. The extraction of the energy spectrum from Eq. (§) is 
routinely performed by Fourier transform, 



<t(w) = I c(t) e iuJt dt. (5) 
Expanding formally ip(t) in terms of eigenf unctions 

m = Y, a nl m §nlme lE " t/h (6) 
nlm 

and inserting Eq. (g) into Eq. (^]) one sees that 

cr(u)=J2^-E n /h)b n . (7) 

n 

Hence, the power spectrum <j(uj) exhibits peaks at the eigenenergies of the system with 
weights given by 
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bn = J2\ a nlm\ 2 = £ | <#)) | <t> n * m ) f (8) 
Im Im 

which are determined by the overlap with the initial wavepacket ip(0). For a finite propaga- 
tion time t the peaks will have a finite width T. Based on an idea by Neuhauser and Wall 
T0|| Mandelshtam and Taylor |TTJ have devised the so called filter-diagonalization method 



as an alternative to the Fourier transform for extracting the energy spectrum from a finite 
time signal c(t). Assuming a form 

c(t)=Y,aje iEjt/h (9) 
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with a,j and Ej being complex the filter-diagonalization allows one to extract E n and b n 
directly from Eq. (0). We will use this stable and accurate method to obtain the spectral 
information from the time signal c(t) which has been calculated semiclassically as follows. 
For the initial wavefunction we have taken a normalized Gaussian wavepacket, ip(0) = 
( 7 2 /7r) 3 / 4 / 7 (r,r ,po) with 

/ 7 (r, r , p ) = exp f-y (r - r ) 2 + ^p Q (r - r )^J (10) 

where the inverse width 7 of the wavepacket and its center (p , r ) in phase space determine 
with which weight the hydrogenic eigenfunctions are covered by ^(O), see Eq. (|]). 

The semiclassical propagator according to Herman and Kluk H is formulated as an 
integral over phase space, 

if 7 (r,r',t) = J J d 3 qd 3 p R 7 (p t ,q t ) exp Q-S (p t ,q t/ 

/ 7 (r, q t , pt) /* (r', q, p) , (11) 

where q 4 = q(i) and p t = p(t) are the phase space values at time t of the trajectory started 
at time t = with (q, p) and propagated under the classical Hamiltonian Eq. (|IJ). The 
action S = J (pq — H) dt accumulated along the trajectory enters Eq. (|1J]) as well as the 
probability density of each trajectory R^(p t ,q t ) which contains all four blocks of the 
monodromy matrix, 




M qq {t) M qp (t) 
M pq {t) M pp {t) 




(12) 



The actual form of the probability density depends on a width parameter 7 which determines 



the admixture of the different blocks Mjj, 



l - \ M qq + M pp - * 7 2 hM qp - JL M pq 



1/2 

(13) 



Although this semiclassical propagator is not uniquely defined through its dependence from 
a suitable chosen parameter 7 it has several important advantages over other forms. Firstly, 
it is globally uniformized since at a caustic _R 7 remains always finite. Secondly, and this 
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is of considerable relevance for practical applications, one does not have to keep track of 
Maslov indices. Instead one has to make i? 7 continuous as the radicant crosses the branch 
cut. Inserting Eq. ( |Tl~D and Eq. fllQ|) into Eq. (0) we obtain a particularly simple form for 
the semiclassical correlation function if the width of the initial Gaussian in ip(0) and of the 
propagator itself in i? 7 are chosen to be the same, 



i(t) = pnh) 3 II rf3qrf3p R ^ Vt ' q *) exp (~h S ( Pt, q *0 9l P ' r °' P °) 9 *i ^ q *' P *' r °' p °) ' 

(14) 



where 



# 7 (q. p ; q'> pO = ex P (~x^ q ~ q '^ 2 ~ ~tf ( p ~ p '^ 2 + in^ p + p '^ q ~ q ')J ( 15 ) 

The integrations over r and r' have been carried out analytically which is possible due to 
our choice of the initial wavepacket as a Gaussian. The remaining integral in Eq. ( IDf) is over 
the entire phase space and in practice c sc (t) is calculated by Monte Carlo integration where 
each randomly chosen phase space point (q, p) represents the initial conditions for a classical 
trajectory. It evolves in time under Hamilton's equations generated by the Hamiltonian of 
Eq. (0) and with the values (q t , p t ) entering Eq. (|i"4|). The number of sampling points 
(trajectories to be run) to achieve convergence depends very much on the initial wavepacket 
ip(0). In general it varies between a couple of thousand and a couple of million trajectories. 

Our first objective is to compare our results with earlier one dimensional calculations 
Although we work in three dimensions we can mimic the one dimensional results to 
some extent by choosing a similar initial wavepacket. 

The result with parameters similar to those from || is shown in Fig. ([J). One sees 
excellent agreement concerning the positions of the peaks with quantum mechanics (crosses) 
and small but noticeable deviations of the weights of the states. This observation, as the 
entire figure, is very similar to the findings of [§] and ||. However, we would like to emphasize 
that our result has been obtained from a 'routinely' applied semiclassical propagator without 
explicit regularization or Langer corrections or any other means implemented to deal with the 
Coulomb singularity. These complication, dealt with in occur only if one uses explicitly 
curved linear coordinates where the problem of the order of operators renders the classical- 
quantum correspondence difficult. This becomes obvious if the semiclassical propagator is 



derived from Feynman's path integral, see e.g., Kleinert's book on path integrals | 12fl . Of 
course, the prize one has to pay in order to avoid these complications in, e.g., a radial 
coordinate, is to work in a higher dimensional (cartesian) space as it has been done here. 

However, even in our approach we should regularize trajectories which hit the Coulomb 
singularity directly (impact parameter zero). Fortunately, these 'head on' trajectories are 
of measure zero among all trajectories contained in the initial conditions and with a Monte 
Carlo method they are hardly ever encountered. Even if such a trajectory is selected by 
chance, one can safely discard its contribution to the propagator. 

The direct semiclassical integration is in principle able to reproduce the spectrum even 
for low excitation as can be seen in Fig. (|2|), and, less surprisingly, for medium excitation 
(Fig. (H|)). However, a systematic trend is apparent from these two spectra: The agreement 
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of the weights is much better to the left of the largest peak than to the right. To understand 
this effect we have plotted in Fig. (||) the average angular momentum fraction 
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contained in the weights b n in addition. One sees that good agreement goes along with a 
large fraction of high angular momentum states in the initial wavepacket and vice versa. 

To support this finding we have prepared a different wavepacket with an additional kick 
(initial momentum) perpendicularly to the axis connecting the center of the wavepacket 
and the Coulomb center. This creates a large fraction of high angular momentum states as 
can be seen in Fig. (|4]). The agreement with the quantum power spectrum is in this case, 
covering the same energy window as [§,0, and Fig. (|I|), much better. Naturally, the one 
dimensional radial calculations of |8||| have only / = states and in Fig. (TJ) the average 
angular momentum £ is also low by construction through the initial state. 

Hence, we can conclude that the power spectrum of hydrogen, including the weights, 
can be reproduced semiclassically. While the semiclassical energies E n are generally in 
good agreement with the quantum eigenvalues the semiclassical weights b n are only accurate 
in the limit of large quantum numbers, i.e. if the initial wave packet contains a large 
fraction of high angular momentum states in each degenerate manifold n. This reflects 
the larger sensitivity of the weights described by off-diagonal matrix elements, compared 
to the (diagonal) energies. Seen in a wider context, our result implies the consequence 
that a one dimensional radial quantum problem is not really one dimensional. Rather, it 
is the limit of angular momentum I = in three, or at least two dimensions. Hence, even 
for large quantum numbers n in the radial problem the semiclassical limit is not reached 
since the angular momentum quantum number is zero. The incomplete semiclassical limit 
causes in the case of the hydrogen problem the remaining discrepancies in the purely radial 
semiclassical spectrum compared with the quantum spectrum. One may also view the failure 
of the one dimensional radial WKB treatment for I = even for large quantum numbers n 
as a consequence of this incomplete semiclassical limit. 

In summary, constructing the time correlation function semiclassically in three cartesian 
dimensions with the help of the Herman-Kluk propagator we have demonstrated that the 
singular Coulomb potential can be treated as any other non-singular interaction without 
any special precautions. Moreover, by virtue of our three-dimensional treatment, we could 
clarify the origin of the discrepancies between the quantum and the semiclassical calculation 
restricted to the radial dynamics only. We hope that this result stimulates future applications 
of semiclassical propagator techniques to Coulomb problems. 
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FIGURES 
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FIG. 1. Semiclassical spectrum (full curve) compared to quantum spectrum (crosses) with 
parameters r = (0, 0, 6000 cm.), Po = (0,0,0) and 7 2 = 2/600 2 . 
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FIG. 2. Semiclassical spectrum (full curve) compared to quantum spectrum (crosses) with 
parameters i-q = (0, 0, 20 a.u.), po = (0,0,0) and = 0.1. 
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3. Lower panel: Semiclassical spectrum (full curve) compared to quantum spectrum 
) with parameters r = (0,0, 2500 a.-u.), Po = (0,0,0) and j 2 = 0.0001; upper panel: 
normalized angular momentum according to Eq. (|l6| ) (circles). 
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FIG. 4. Lower panel: Semiclassical spectrum (full curve) compared to quantum spectrum 
(crosses) with parameters ro = (0,0, 6000a.it.), po = (0,0.0075,0) and r y 2 = 2/600 2 ; upper panel: 
average normalized angular momentum according to Eq. fllq) (circles). 
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